Island morphology and adatom self-diffusion on Pt(lll) 



OS 

o> 



Ghyslain BoisvertQ and Laurent J. LewisB 

Departement de Physique et Groupe de Recherche en Physique et Technologie des Couches Minces (GCM), Universite de 
Montreal, Case Postale 6128, Succursale Centre-Ville, Montreal, Quebec, Canada H3C 3J7 

Matthias Scheffler 

Fritz- Haber-Institut der Max- Planck- Ges ells chaft, Faradayweg 4-6, D-1^195 Berlin- Dahlem, Germany 

Submitted to Physical Review B 

February 1, 2008 



(N 



G 

O 

o 



> 

(N 

in 
o 

OS 
-i— > 



i 

O 
o 



X 



The results of a density-functional-theory study of the formation energies of (100)- and (lll)-faceted 
steps on the Pt(lll) surface, as well as of the barrier for diffusion of an adatom on the flat surface, 
are presented. The step formation energies are found to be in a ratio of 0.88 in favour of the 
(lll)-faceted step, in excellent agreement with experiment; the equilibrium shape of islands should 
therefore clearly be non-hexagonal. The origin of the difference between the two steps is discussed 
in terms of the release of stress at the surface through relaxation. For the diffusion barrier, we also 
find relaxation to be important, leading to a 20% decrease of its energy. The value we obtain, 0.33 
eV, however remains higher than available experimental data; possible reasons for this discrepancy 
are discussed. We find the ratio of step formation energies and the diffusion barrier to be the same 
whether using the local-density approximation or the generalized-gradient approximation for the 
exchange-and-correlation energy. 
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I. INTRODUCTION 

Detailed knowledge of surface properties is important 
to understanding a wide variety of phenomena such as 
catalysis, surface reactivity, growth, etc. Of particular 
importance are such quantities as the step formation en- 
ergy (SFE), which determines the equilibrium shape of 
islands on flat terraces, and the potential energy surface 
seen by, e.g., an adatom, which provides information on 
the preferred sites for adsorption and the kinetics of diffu- 
sion. We present here a detailed study of these properties 
for the (111) surface of platinum within the framework 
of density- functional theory (DFT).EJ 

The (111) surface of fee metals, and in particular Pt, 
is of interest for (at least) two reasons. First, as dis- 
played in Fig. [l], it exhibits two possible step geome- 
tries, named according to the micro-facet that step-edge 
atoms form with atoms in the layer underneath, namely 
a (lOO)-faceted step, where an edge atom has a single 
near neighbor at the base of the step (i.e., the micro- 
facet constitutes a square lattice) and a (lll)-faceted 
step, where each atom has two neighbors (triangular lat- 
tice). The only difference between the two geometries, 
as far as nearest-neighbors are concerned, is that atoms 
at the base of the step have a coordination of 10 for 
the (lOO)-faceted step and 11 for the (lll)-faceted step. 
The two steps are evidently very similar and the forma- 
tion energies are thus expected to be comparable, i.e., the 
equilibrium island shape should be very nearly hexago- 



nal, with the sides consisting of, alternately, (100)- and 
(lll)-facgted steps- This has indeed been observed in the 
case of Iid and Ag.B For Pt, however, scanning-tunneling 
microscopy (STM) reveals, rather, a strong preference 
for (lll)-faceted steps — 0.87±0.02 as measured by the 
ratio of step formation energies per unit length.a On 
the theory side, DFT calculations in the local-density 
approximations (LDA) by Stumpf and Scheffler predict 
this behaviour in the case of Al, with a ratio of SFE 
of 0.93;Eil this has not yet been confirmed experimentally. 
For Pt, in contrast, corresponding calculations by Feibel- 
man failed to reproduce the experimental results, leading 
to essentially equal SFE for the two types of steps .0 

A second feature that makes the fcc(lll) surface inter- 
esting is that it possesses two different adsorption sites: 
the fee (normal) site, where an adatom sits in a posi- 
tion appropriate to the stacking of atomic planes in a 
perfect fee crystal, and the hep (stacking fault) site, cor- 
responding to an hep stacking. Both sites have three 
nearest neighbors, the difference between the two lying 
in the second-layer arrangement, as illustrated in Fig. |^. 
It has been shown that, for transition metals, the pre- 
ferrecUsite for adsorption depends on the filling of the d 
bandMEl For Pt, DFT-LDA calcaila.tions predict the fee 
site to be much more favourable, Q't2l by 0.12-0.18 eV, de- 
pending on the state of relaxation of the substrate. This 
is the largest difference obseispdi-^ o far) for late transi- 
tion metals and noble metals .Br HI 12 ! Experimentally, also, 
there is evidence that the fee site is preferred over the 
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hep site in Pt;t2H13 while a precise numerical value can- 
not be inferred from the measurements, it is-,estimated 
that the difference should be at least 0.06 eVila However, 
despite this,japparent agreement, the DFT calculations 
for Pt(lll)mL9 have been unable to reproduce correctly 
the diffusioit-hai'rier — 0.38-0.41 versus ~0.26 eV from 
experiment JMj'EJ Since it is of primary importance to have 
reliable and accurate energy barriers in order to predict 
growth (see for example Ref. [l5|) , it is essential that this 
discrepancy be resolved. 

Clearly, a quantitative picture of the surface proper- 
ties of Pt(lll) is still missing. In order to address this 
problem, we present here the results of extensive ab-initio 
total-energy calculations of the formation energies of the 
two kinds of steps, as well as of the energetics of adatom 
adsorption and diffusion. The calculations have been car- 
ried out within the LDA, but we have also carried out a 
series of calculations in the generalized-gradient approx- 
imation (GGAjlj so as to assess the applicability of the 
LDA to this system. We find the LDA to offer a bet- 
ter description (compared to experiment) of the lattice 
constant of bulk Pt tha& the GGA, in agreement with 
previous calculationsJlU'EJ as well as some properties of 
the clean (111) surface (surface relaxation and work func- 
tion). Within numerical accuracy, however, we observe 
no sizeable effect on the ratio of SFE and on the energet- 
ics of adatom adsorption and diffusion. 

As discussed in more detail below, the electronic wave- 
functions were expanded in plane waves, which enables 
us to deal easily, and completely, with the effect of relax- 
ation on adsorption and adatom diffusion. In the calcula- 
tions of Refs. and [l^, only nearest-neighbor relaxation 
was included, at best. We observe that proper account 
of atomic relaxation leads to a significant decrease, by 
^20%, of the value of the diffusion barrier on the flat 
(111) surface. .Qui. estimate for this quantity, while. closer 
to experimentEJo than previous calculations,HL9 how- 
ever remains high — 0.33 versus ~ 0.26 eV; possible rea- 
sons for this discrepancy will be discussed. In contrast, 
the difference in energy between the two adsorption sites 
increases upon relaxing, from 0.10 to 0.17 eV. 

One other important advantage of plane waves is that, 
contrary to the Gaussian orbitals used in Ref. 0, they are 
independent of the atomic positions and should there- 
fore provide a more adequate description of the subtle 
difference in SFE expected here. In agreement with ex- 
periment, our highly-converged calculations indicate a 
clear preference for (lll)-faceted steps over (lOO)-faceted 
steps, the SFE being, in a ratio of 0.88 (versus about 
0.87 experimentally) .□ Our calculations, further, provide 
a simple explanation for the origin of the energy differ- 
ence between the two steps in terms of the release of 
surface stress through relaxation. Before discussing our 
results in detail, we give a brief description of our com- 
putational approach. 



II. COMPUTATIONAL DETAILS 

As already noted above, the calculations reported 
here were performed within the framework of density 
functional theoryjd using both the LDAEl and the GGAE3 
for the exchange-and-correlation energy. The ion cores 
were approximated by pseudopotentials with 5d electrons 
treated as valence states. The pseudopotentials were gen- 
erated using the semi-relativistic scheme of Troullier and 
Martinaij and expressed in the KleinmanpBy-lander form 
using the s component as the local one.E3~E3 The elec- 
tronic wave-functions were represented using a plane- 
wave basis set with kinetic energy up to 40 Ry in the 
LDA and 45 Ry in the GGA. To improve convergence, 
the electronic states were occupied according to a Fermi 
distribution with ksT c \ = 0.1 eV and the total energies 
obtained by extrapolating to zero electronic temperature. 
For similar reasons, the calculations were initiated using 
wave-functions obtained from the self-consistent solution 
of the Kohn-Sham Hamiltonian in a mixed basis set com- 
posed af-pseudo-atomic orbitals and plane waves cut off 
at 4 RyE3 The minimization of the energy with respect to 
the electronic degrees of freedom was done using an itera- 
tive procedure.EJ After achieving electronic convergence, 
the atoms were moved according to a damped Newton 
dynamics until forces became less than 0.01 eV/A. All 
the calculations were performed using the supercell ap- 
proach. Details of the cell shape and size, as well as k- 
point sampling, for the different geometries considered, 
are given along with the results in the following sections. 

III. RESULTS 

A. Bulk and Surface Properties 

In order to assess the validity of our approach, and for 
completeness, we first determined the lattice constant of 
the bulk material as well as the properties of the clean 
Pt(lll) surface — surface energy, relaxation, excess sur- 
face stress, and work function. As a reminder, the surface 
stress tensor, g a p, is given by 

g a (3 = 7<W + d~f/de a/ 3 (1) 

where 7 is the surface energy per unit area, £ a p is the 
strain tensor and 5 a p is the Kronecker delta function. 
The second term in this equation represents the excess 
surface stress. For the clean surface, we used a (1 x 1) 
supercell consisting of 5 or 7 (111) atomic planes plus 
~io A of vacuum. Integration of the first Brillouin zone 
was done over a uniform grid of 100 k points in the x — y 
plane, which was found to yield well-converged results, 
e.g., within 2 meV for the surface energy. The latter was 
calculated by comparing to a bulk-like (lxl) supercell 
containing 3 layers (and of course no vacuum). The same 
k-point density was used for the x — y plane; to compen- 
sate for the smaller size of the cell in the z direction, 
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the two-dimensional grid was replicated 4 times along z 
so as to get a density of points similar to that in the 
x — y plane. For the excess surface stress, we varied the 
in-plane lattice constant for both the clean surface and 
the bulk, while keeping the atoms' z coordinates fixed 
to their bulk-like values, and examined the concomitant 
variations in the total energy. 

The results are listed in Table |, along with those from 
other ab-initio calculations and experimental values when 
available. Evidently, our LDA lattice constant is con- 
sistent with previous calculations and with experiment. 
The GGA value, in contrast, while in agreement with 
other calculations, overestimates somewhat the lattice 
constant, by more than 2%; it is a well-known fact that 
the GGA yields larger lattice constant than the LDAta 
— actually overcompensates in the present case. For the 
surface energy and excess surface stress, our LDA results 
agree well with previous calculations. We note that the 
excess surface stress is large and positive, meaning that 
the Pt(lll) surface is under significant tensile stress, i.e., 
would prefer a smaller lattice constant. Unfortunately, to 
our knowledge, there exists no experimental determina- 
tion of these quantities. In the GGA, both the surface 
energy and the excess surface stress decrease. This effect 
of the GGA on the surface-energy was actually predicted 
from jellium calculations!!^ 

We have also calculated the top and second layer re- 
laxation, Ae?i2 and Ac?23, i-e., the change in interlayer 
spacing relative to the bulk value. We find a small (0.4%) 
outward relaxation for the top layer. This is quite a bit 
smaller than the value reported in Ref. [lO], 1.25%. The 
origin of the discrepancy between the two LDA calcula- 
tions is not clear; it may be due to different choices of 
basis sets (LCAO in Ref. [lO] versus plane waves here). 
The GGA, interestingly, leads to a small inward relax- 
ation. Currently available experimental data vary widely, 
in the range 0-2.5%, and are therefore not of much help in 
resolving the issue. For the second layer, theory and ex- 
periment agree that it should be insignificant, i.e., Ag?23 
is a small fraction of a percent. 

For the work function, finally, the experimental values 
also vary quite a bitpip the range 5.77-6.10 eV. According 
to Kaack and Fickjlj however, the work function has a 
small temperature dependence, decreasing slightly with 
temperature. Since our calculations are performed at 
K, we expect that they should compare well with the 
largest experimental values. We find, indeed, that the 
LDA result, 6.07 eV, is in excellent agreement with the 
largest experimental number, 6.10 eV. The GGA value, in 
contrast, is significantly smaller — 5.70 eV — indicating, 
once more, that the LDA provides a better description of 
Pt than the LDA. In spite of this, the two approximations 
will carefully be examined in the context of SFE and 
diffusion barriers. 



B. Step Formation Energy 

We come now to the heart of the matter, namely the 
energetics of step formation. In order to determine the 
SFE, wc constructed vicinal surfaces (using rectangular 
surface cells) appropriate to each type of steps. For the 
(lOO)-faceted step, we examined both a (211) and a (332) 
surface; the former has 3 atoms per terrace while the 
latter has 5. For the (lll)-faceted step, only the (221) 
surface, which contains 4 atoms per terrace, is considered. 
Again, here, a vacuum region of approximately 10 A was 
included in all cases. 

The energies of the vicinal surfaces were determined 
using the same procedure as in the case of the clean sur- 
face, i.e., by comparing to an appropriate bulk model. In 
order to minimize the error arising from the use of differ- 
ent geometries, the bulk reference system for a given vic- 
inal surface was always taken to have the same in-plane 
geometry as the surface. Thus, the same k points were 
used in the x — y plane, while for the z coordinate, the 
grid was adjusted to yield a comparable density. In to- 
tal, the bulk supercells corresponding to the (211), (322), 
and (221) surfaces contained 6, 34, and 18 atoms, respec- 
tively. 

The SFE is given, simply, by the difference in energy 
between a surface with a step and one without. Since we 
are dealing with vicinal surfaces here, this is equivalent 
to subtracting from the vicinal-surface energy (per ter- 
race), ovic, that portion of the clean-surface energy cor- 
responding to the exposed (111) area.Q If we call C(m) 
the clean-(lll)-surface energy per atom and neglect step- 
step interactions we find, for the (lOO)-faceted step: 



^(100) 



(N 



and for the (lll)-faceted step: 



-fr(lll) 



(N 



(2) 



(3) 



where N is the number of atoms per terrace, as defined 
earlier. 

It is evident from Eqs. [2] and || that, in order to de- 
termine reliably the SFE ratio, very accurate surface en- 
ergies are required for both vicinal and clean surfaces. 
It is our purpose here to assess carefully the accuracy 
of our calculations through a detailed convergence study. 
As explained earlier, the error arising from the supercell 
geometry is minimized by always comparing surface and 
bulk energies obtained using the same in-plane period- 
icity and k-point density. Of course, it is essential that 
the energies be converged with respect to Brillouin-zone 
integration; for the clean (111) surface, the k-point sam- 
pling scheme used here leads to values converged within 
2 meV, as discussed in Sec. |III A . For the vicinal sur- 
faces, we used a similar sampling scheme and, as we will 
see below, the error is of the order of a few meV, so that 
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differences in energy (e.g., between steps) of a few hun- 
dredths of an eV are significant. 

The results for the (lOO)-faceted step under a variety 
of theoretical conditions are listed in Table O. First, we 
examine the effect of relaxation and size within the LDA. 
With all atoms in the bulk-like configuration (referred to 
as "rigid" in Table [FJ), we find the SFE to be rather in- 
sensitive to the size of the supercell: adding two layers to 
the (211) slab increases the SFE by a mere 0.02 eV/ (step 
atom), while using a (322) surface instead of a (211) leads 
to a small decrease of 0.01 eV/ (step atom). (The number 
of layers refers to the number of (lll)-like layers in the 
slab before a rotation is applied to make the surface vic- 
inal.) However, upon relaxing all atoms, except those in 
the central (111) layer (configurations referred to as "re- 
laxed" in Table 0) , the SFE is found to decrease strongly, 
from 0.62 to 0.43 eV/(step atom) for the (322) surface. 
For the relaxed configurations, we have also examined the 
convergence with respect to the k-point sampling. In all 
cases, the SFE changes by at most 0.01 eV/(step atom) 
upon increasing the number of k points. Thus, we esti- 
mate the SFE for the (lOO)-faceted step to be 0.43±0.02 
eV/(step atom) within the LDA. 

For the (lll)-faceted step, in view of the above results, 
we have studied a single vicinal surface, namely the (221), 
at fixed and converged k-point density, as indicated in 
Table [II. Again, here, relaxation affects strongly the 
SFE, which decreases very markedly — from 0.64 to 0.38 
eV for a 5-layer slab in the LDA. However, increasing the 
thickness from 5 to 7 layers brings about no significant 
changes in the SFE. Thus, our best LDA-SFE value for 
the (lll)-faceted step is 0.38±0.02 eV/(step atom). 

Our calculations indicate, therefore, that the (111)- 
facetcd step has a lower formation energy than the (100)- 
facctcd step 0.38 versus 0.43 eV (in the LDA), leading 
to a ratio Ef^JE?*^ of 0.88±0.07. This is in excellent 

agreement with the experimental value of 0.87±0.02,§ 
but at variance with a previous LDA calculation by 
Feibelmanj3 who found that the two steps are nearly 
equivalent, that is, 0.46 eV for the (lll)-faceted step 
versus 0.47 eV for the (100), i.e., a ratio of 0.98. The 
SFE values differ from Feibelman's not only in a rela- 
tive sense, but also in an absolute sense: the values we 
find are significantly smaller, by 0.08 eV/(step atom) for 
the (lll)-faceted step and 0.04 eV/(step atom) for the 
(lOO)-faceted step. Though the reasons for these differ- 
ences are not clear, they may originate in the choice of 
basis sets: while we use plane waves, Feibelman employs 
Gaussian orbitals which are more sensitive to the details 
of the atomic configuration, as discussed in Ref. 0. In 
view of this, it might perhaps be the case that a Gaus- 
sian basis set lacks the accuracy needed to resolve such 
small energy differences as those involved here. 

A value of 0.37 eV/(step atom) for the (lll)-faceted 
SFE has also been estimated from the experimental sur- 
face freejienergy anisotropy between the (110) and (111) 
surfaces.cd While this corresponds quite closely to our 



value of 0.38 eV/(step atom), the agreement is fortuitous 
since the above result was obtained assuming a value of 
the surface energy of 0.097 eV/A 2 , much lower than that 
calculated here, 0.124 eV/A 2 . [Using the latter value for 
the surface energy would lead to a (lll)-faceted SFE 0.45 
eV/ (step atom) in the approach of Ref. |d].] 

Recently, some concerns have been expressed regard- 
ing the procedure used here to determine the surface en- 
ergy, whiph | _should diverge as the thickness of the slab 
increases OO We reinterpreted our results using the ap- 
proach suggested in Ref. ^3| and found only small changes 
in the (HI) surface energy, now 0.121 eV/A 2 rather than 
0.124 eV/A 2 ; for the SFE, we obtain now 0.41 and 0.46 
eV/(step atom) for the (HI)- and (lOO)-faceted step, re- 
spectively, compared to 0.38 and 0.43 eV/(step atom) 
using the usual approach. The SFE ratio remains ap- 
proximately unchanged, 0.89 versus 0.88. We are thus 
led to conclude that, while the uncertainty on the SFE 
might be of the order a 0.03 eV, the value we find for the 
ratio is accurate to a few percent, and is not affected by 
the numerical procedure used. 

As mentioned in the previous section, the LDA seems 
to provide a better description of bulk Pt, as well as of 
the (111) surface, than the GGA. In view of this, it is 
expected that the vicinal surfaces are also better rep- 
resented within the LDA. The question remains open, 
however, because there exists no firm experimental data 
to compare our results to, and it is therefore of interest 
to calculate the SFE also within the GGA. The results 
are given in Tables H and |l| We observe the GGA-SFE 
to be systematically lower than the corresponding LDA 
values, as is also true of the (111) surface energy. Fur- 
ther, convergence with respect to both size and k-point 
density is similar in the two approximations. We there- 
fore conclude to GGA values of 0.29±0.02 and 0.25±0.02 
eV/(step atom) for the (100)- and (lll)-faceted step, re- 
spectively. The resulting SFE ratio is 0.86±0.10, essen- 
tially unchanged from the LDA value, namely 0.88±0.07. 

In order to understand why the two steps have different 
formation energies, it is of interest to consider, first, the 
(111)- to (lOO)-faceted SFE ratio in the unrelaxed (bulk- 
like) configuration. We find, from Tables [0] and III, this 
ratio to be equal to 1.03±0.07 (using the LDA), compared 
to about 0.88 for the relaxed configurations, as we have 
seen above. Thus, before the atoms relax, the two steps 
are nearly equivalent [with perhaps a slight preference for 
the (lOO)-facctcd step], as could be expected from a sim- 
ple nearest-neighbor model as explained in the Introduc- 
tion. Evidently, therefore, the observed step anisotropy is 
closely related to relaxation, and this can be u nders tood 
in the following way: As we have seen in Sec. Ill A, the 
Pt(lll) surface is under large tensile stress, which can 
be locally relieved at steps. However, because the atomic 
configurations are different (albeit slightly) for the two 
kinds of steps, the relaxation patterns also differ, and 
lead to different energetics. This can in fact be seen very 
clearly in Fig. ^ where we plot the displacement pat- 
terns for the two types of steps: Some atoms suffer very 
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large displacements — by as much as a few percent (rela- 
tive to the bulk nearest-neighbor distance) for those that 
sit closest to the steps. More important, it is also clear 
from this figure that the displacements associated to the 
(lll)-faceted step are larger than for the (lOO)-faceted 
step, i.e., the former can relieve stress more efficiently 
than the latter, and is thus energetically more favorable. 
It should be mentioned that the displacements we find 
here differ from Feibelman'sQ by as much as 1% in some 
cases, and might possibly explained the discrepancy be- 
tween the two sets of results; this is likely related, again, 
to different choices of basis functions. 

The relation between relaxation and stress can be un- 
derstood in a more quantitative manner by considering 
the change in energy resulting from the displacement in- 
wards of an edge atom, i.e., in the direction normal to the 
step and parallel to the terrace. Starting with both step 
models in their bulk-like geometry and moving an edge 
atom by the same amount for the two steps, we find that 
the (lll)-faceted step, because of its triangular geometry, 
releases more energy than the (lOO)-faceted step: for a 
displacement of 0.14 A, corresponding approximately to 
the observed relaxation, we find the (lll)-faceted step to 
be already 22 meV lower in energy than the (lOO)-faceted 
step. 

Clearly therefore, the difference in SFE arises from the 
large excess surface stress of the Pt(lll) surface — 0.25 
eV/A 2 (cf. Table |h, about twice as large as the surface 
energy, 0.124 eV/A 2 . We may compare this with the 
corresponding situation for Ir/Ir(lll), where the equi- 
librium island shape is nearly hexagonaLM. In this case, 
the excess surface stress is 0.128 eV/A 2 ] 2 ^ significantly 
smaller than the surface energy, 0.204 eV/A 2 . Evidently, 
large surface energies lead to large SFE, and large excess 
surface stresses to large differences between the two steps. 
These observations therefore suggest that, as a "rule of 
thumb", the SFE ratio should differ from one (i.e., non- 
hexagonal equilibrium island shape) when the excess sur- 
face stress is sizeably larger than the surface energy. This 
is in fact the case of Al/llip ami Au(lll),E3 while the 
opposite is true of Rio and Cuo To our knowledge, no 
information on the equilibrium island shape is available 
for Rh and Cu, but we would predict it to, be hexago- 
nal in both cases. For Au, reconstructioncH is likely to 
be important in determining the island shape. For Al, 
finally, ab initio calculations of the kind presented here 
have been performed by Stumpf and Schefflerja and a 
SFE ratio of 0.93 is indeed found. In this case, the ra- 
tio does not seem to be affected by relaxation, contrary 
to our results for Pt, but different electronic orbitals are 
involved — sp for Al and d for Pt. 

There have been other calculations of the equilibrium 
island shape on Pt(lll). Using a tight-binding model, 
Papadia et al. found essentially no difference between 
the two steps in their bulk-like configuration—as is the 
case here, but did not consider relaxationo Within 
equivalent-crystal theory Khare and Einstein found a ra- 
tio close to unity (0.968) for the SFE, without allowing 



in-plane relaxation;^ also, in this approach, the surface 
energy is predicted to be 0.076 eV/A 2 , quite a bit smaller 
than our 0.124 eV/A 2 . Fully- relaxed calculations using 
the semi-empirical embedded-atom method (EAM) have 
also been performed and lead to a SFE ratio very close 
to one;E3 it is however doubtful that the EAM potential 
is robust enough to account for the small energy differ- 
ences involved here. Other approaches have been pro- 
posed, -based on coordination- or orientation-dependent 
bondSjOL 2 ] which do not take relaxation effects into ac- 
count. 

We also list, in Tables || and III, the work function 
for the different surfaces examined. As expected (see fog 
instance Ref. |33|), and already observed by Feibelman,Q 
the work function is smaller for the vicinal surfaces than 
for the (111) surface. We also observe, in the case of the 
(lOO)-faceted step, a small dependence on the the terrace 
length. This agrees with Feibelman's calculations, whiki 
a stronger dependence is reported from experiment .lj 
This might be due to the fact that terraces studied in 
experiment are. much wider than ours and/or surface 
contamination.0 



C. Atom Adsorption and Diffusion 

We now discuss adsorption and diffusion of a Pt atom 
on the Pt(lll) surface. For these calculations, an adatom 
is added on one surface of the Pt slab while the other 
surface is constrained to its bulk-like configuration. In 
order to determine the energies at the two adsorption 
sites as well as the barrier for diffusion, we considered 
both a (2 x 2) and a (3 x 3) cell with, again, approximately 
10 A of vacuum. Unless otherwise noted, the integration 
over reciprocal space was performed using a mesh of 16 
equidistant k points for the (2 x 2) cell and 9 k points 
for the (3 x 3) cell. 

First, starting with the (2 x 2) system, we examined 
convergence with respect to the number of layers, which 
we varied from 3 to 6. The results are given in Table 
IV . We observe significant changes upon going from 3 to 
4 (for a given state of relaxation), while increasing this 
number further brings about changes of at most 0.01 eV; 
thus, 4 layers seem to be sufficient for reliable estimates 
of both the diffusion barrier, Ed, and the difference in 
adsorption energies, AE a d s - 

The barrier we obtain for the unrelaxed substrate, 0.41 
eV, and the difference in adsorption energies, 0.10 eV, 
agree well with previous LDA calculations.Qil3 It is how- 
ever clear from Table IV that relaxation effects are again 
here important: If we allow the topmost layer to relax, 
Ed decreases to 0.34 eV and Ai? a( js increases to 0.17 eV; 
including second-layer relaxation as well results in rela- 
tively minor changes to the energies (less than 0.02 eV). 

It should be noted that, because of the asymmetry 
between the two adsorption sites, there are in fact two 
barriers for diffusion. However, the difference between 
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the two barriers is such that, for temperatures of inter- 
est, it is the highest-energy barrier that limits diffusion, 
i.e., the adatom will escape rapidly from the low-energy 
adsorption state but get trapped in the high-energy site. 

We have also examined convergence with respect to 
k-point sampling, en ergy cutoff, and lateral size. As in- 
dicated in Table IV, we find in all cases very modest 
changes of at most 0.01 eV. Likewise, using the GGA 
does not lead to appreciable changes to the energy bar- 
rier, while Ai? a ds decrea ses b y about 0.03 eV. For the 
reasons discussed in Sec. [II A, we suspect that the LDA 



values are more accurate; our best, highly-converged esti- 
mates of E d and AE ads are thus 0.33±0.03 and 0.17±0.03 
eV, respectively. As discussed in Ref. g, the large differ- 
ence between the fee and the hep site is related to the an- 
gular character of the d orbitals, whose bonding strength 
depends on the filling and radial quantum number of the 
bands. Due to this difference, the transition site for jump 
diffusion does not lie exactly midway between the two 
equilibrium sites but, rather, about 0.07 A towards the 
hep site. The potential energy surface is however very 
flat in this region, changing by no more than 0.01 eV 
upon going from the transition to the midpoint site. 

There is very little experimental information available 
for Ai5 a d s ; as mentioned in the Introduction, only a lower 
limit of 0.06 eV has been determined,E2 and this is consis- 
tent with our results. For the diffusion barrier, a value of 
0.25±0.02 eV has been inferred from field-ion microscopy 
(FIM) measurements in the temperature range 92-100 
K.Ej Also, based on a comparison between STM mea- 
surements of the island density and kinetic Monte Carlo 
simulations between 110 K, and 160 K, an estimate of 
0.26±0.01 eV is obtained.ta Both values agree, but dis- 
agree with our result of 0.33±0.03 eV, which is surprising 
in view of the high-level of accuracy and convergence of 
our calculations. 

One possible explanation for this disagreement would 
be the neglect, in our calculations, of dynamical effects. 
These cannot be assessed directly from first-principles 
but have been shown to be insignificant at low temper- 
atures, either from the calculation pf dynamical correc- 
tions to the transition state theptvEj or from empirical 
molecular-dynamics simulations. LjEj Experimentally, it 
is clear that the FIM valuald suffers from poor statistics 
— only two points are used to determine the Arrhenius 
parameters — and therefore the error bar is large. Con- 
cerning the STM experiment, it has been pointed out by 
the authors that small-cluster mobility could affect the 
Monte-Carlo estimate of the barrier if its energy is close 
to that for adatom diffusion, in which case the quoted 
value would be a lower bound to the actual barrier. As 
a final point, it should be mentioned that many different 
empirical potentials have been used to delepnpe the dif- 
fusion barrier of Pt adatoms on Pt(lll) Jl3l£j l£3 yielding 
to values in the range 0.01-0.18 eV, i.e., much lower than 
the experimental value, which we argue, is a lower bound 
to the actual barrier. Thus, such models are clearly too 
crude to provide a proper description of the energetics of 



diffusion for the present system. 



IV. SUMMARY 

We have used highly-accurate ab-initio methods to cal- 
culate the ratio of (111)- to (lOO)-faceted step formation 
energies on the (111) surface of Pt. We find, in excel- 
lent agreement with experiment, (lll)-faceted steps to 
be favoured over (100) in a ratio of about 0.88; as a 
consequence, the island on this surface should be clearly 
non-hexagonal. The difference between the two steps is 
related to the large tensile stress of the Pt(lll) surface, 
which is released in a different manner because of differ- 
ences in the local topology. Our calculations underline 
the importance of relaxation in such cases: while the 
two steps are about equivalent for the unrelaxed sub- 
strate, relaxation does bring about large changes in the 
formation energies. Likewise, relaxation is important to a 
proper determination of equilibrium-site energies on the 
flat (111) surface. We find the fee site to be preferred 
over the hep by a sizeable 0.17 eV (after full relaxation), 
consistent with experiment which provides a lower bound 
of 0.06 eV for the difference between the two sites. 

We have also calculated the energy barrier for adatom 
diffusion and found a fully-relaxed value of 0.33 eV. While 
constituting an improvement over previous calculations, 
this value remains larger than experiment, by about 0.07 
eV. The discrepancy might be due to limitations of our 
theoretical approach (e.g., finite size), but it might also 
due to errors in the interpretation of the experimental 
data — poor statistics, neglect of small-cluster contri- 
butions (e.g., dimers) to mass transport, i.e., incorrect 
assumption regarding the critical nucleus size. More ex- 
periments are needed to clarify this point. Likewise, it 
would be of interest that measurements of the SFE and 
difference in energy between fee and hep sites be carried 
out so as to assess the validity of LDA (versus GGA) in 
the present context. 
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TABLE I. Bulk and clean (111) surface properties of platinum: lattice constant a, surface energy a, excess surface stress t, 
top and second layer relaxation, Adi2 and Ad23, and work function W . 





a 




T 


Ad 12 


Ad 23 


w 




(A) 


(eV/A 2 ) 


(eV/A 2 ) 


(%dbulk) 


(%dbulk) 


(eV) 


Present LDA 


3.92 


0.124 


0.25 


0.4 


-0.2 


6.07 



Present GGA 
Other LDA 3.87 
Other GGA 
expt. 




0.213 



If 0.2890 



<| 2.5 |P'B< 
1.1±0.5E3 



l.oB < 



0.0 r 



<| 0.4 [P 



5.7B 
6.10E 



-0.050 
0.00 5.77-6.1 



TABLE II. Step formation energy (SFE) for the (lOO)-faceted step on Pt(lll) (in eV per step atom) and work function W 



for the corresponding vicinal surfaces (see text), as a function of the number of (111) layers, Nl, number of k points, Nk, and 
approximation scheme for the exchange-and-correlation energy (XC). 



surface 


N L 


N k 


XC 


SFE 


W 










(eV/at) 


(eV) 


(211)-rigid 


5 


8 


LDA 


0.63 


5.88 


(211)-relaxed 


5 


8 


LDA 


0.45 


5.92 


(211)-relaxed 


5 


12 


LDA 


0.46 


5.93 


(211)-rigid 


7 


8 


LDA 


0.65 


5.89 


(211)-relaxed 


7 


8 


LDA 


0.47 


5.94 


(211)-rigid 


5 


8 


GGA 


0.50 


5.58 


(211)-relaxed 


5 


8 


GGA 


0.34 


5.63 


(211)-relaxed 


5 


12 


GGA 


0.34 


5.64 


(211)-rigid 


7 


8 


GGA 


0.52 


5.59 


(211)-relaxed 


7 


8 


GGA 


0.35 


5.63 


(322)-relaxed 


5 


4 


LDA 


0.44 


5.96 


(322)-rigid 


5 


8 


LDA 


0.62 


5.94 


(322)-relaxed 


5 


8 


LDA 


0.43 


5.95 


(322)-rigid 


5 


4 


GGA 


0.46 


5.64 


(322)-relaxed 


5 


4 


GGA 


0.29 


5.66 



TABLE III. Same as Table II but for the (lll)-faceted step. 



surface 


N L 


N k 


XC 


SFE 
(eV/at) 


W 

(eV) 


(221)-rigid 


5 


8 


LDA 


0.64 


5.84 


(221)-relaxed 


5 


8 


LDA 


0.38 


5.80 


(221)-relaxed 


7 


8 


LDA 


0.38 


5.82 


(221)-rigid 


5 


8 


GGA 


0.50 


5.54 


(221)-relaxed 


5 


8 


GGA 


0.25 


5.50 







TABLE IV. Adatom diffusion barrier Ed and difference in adsorption energies between fee and hep sites, 
A£ ads = E^ s - E^, for Pt on Pt(lll) under various calculational conditions, as discussed in the text. 



xc 


supercell 


allowed to relax 


k point grid 


E d 


A-E a ds 










(eV) 


(eV) 


LDA 


(2x2), 3- layer 


adatom only 


4x4 


0.47 


0.22 


LDA 


(2x2), 3- layer 


adatom + top layer 


4x4 


0.41 


0.26 


LDA 


(2x2), 4- layer 


adatom only 


4x4 


0.41 


0.10 


LDA 


(2x2), 4- layer 


adatom + top layer 


4x4 


0.34 


0.17 


GGA 


(2x2), 4- layer 


adatom + top layer 


4x4 


0.33 


0.14 


LDA 


(2x2), 4- layer 


adatom + 2 top layers 


4x4 


0.36 


0.17 


LDA 


(2x2), 4- layer 


adatom + top layer 


5x5 


0.35 


0.17 


LDA(50 Ry) 


(2x2), 4- layer 


adatom + top layer 


4x4 


0.35 




LDA 


(2x2), 5- layer 


adatom + top layer 


4x4 


0.33 




LDA 


(2x2), 6- layer 


adatom + top layer 


4x4 


0.35 




LDA 


(3x3), 4- layer 


adatom + top layer 


3x3 


0.33 
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(a) (100) -faceted step 



(111) -faceted 




FIG. 1. Top view of an island (filled circles) on a fcc(lll) 
surface indicating the two possible step geometries. 



top 



view 



side 
view 
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FIG. 2. Top and side view of the two different adsorption 
sites on a fee (111) surface: fee site (filled circle) and hep 
site (dashed circle). The side view corresponds to the plane 
indicated by the straight line in the top view. 



(b) (111) -faceted step 




FIG. 3. Side view of the supercell showing the relaxed 
geometry for (a) the (lOO)-faceted step [as obtained from the 
(322) surface] and (b) the (lll)-faceted step [from the (221) 
surface]. The filled symbols represent the middle (111) layer, 
which is fixed in its bulk-like position. The circles are in the 
plane of the sheet while the squares are in a different plane, 
distant by d/2, where d is the nearest-neighbor distance. The 
arrows indicate the periodicity of the supercell. 
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